arXiv: 1501.06737vl [astro-ph.HE] 27 Jan 2015 


On the origin and physics of Gamma Flares in Crab Nebula 


George Machabeli, Andria, Rogava and David Sha.pakidze 
Centre of Theoretical Astrophysics, Institute of Theoretical Physics, Ilia State University, 

Tbilisi, Georgia 


Received 


accepted 


Submitted to Astrophysical Journal 




ABSTRACT 


We consider parametric generation of electrostatic waves in the magneto¬ 
sphere of the pulsar PSR0531. It is shown that in the framework of this mech¬ 
anism it is possible to convert the pulsar rotational energy into the energy of 
Langmuir waves. The maximum growth rate is achieved in the “superluminal” 
area, where phase velocity of perturbations is exceeding the speed of light. There¬ 
fore electromagnetic waves do not damp on particles. Instead, they create plas- 
mon condensate, which is carried out, outside of the pulsar magnetosphere and 
reaches the Crab nebula. It is shown, that the transfer of the energy of the plas- 
mon condensate from the light cylinder to the active region of the nebula happens 
practically without losses. Unlike the plasma of the magnetosphere, the one of 
nebula contains ions, i.e., it may sustain modulation instability, which leads to 
the collapse of the Langmuir condensate. Langmuir wave collapse, in turn, leads 
to the acceleration of the distribution function particles. Furthermore, we con¬ 
sider processes leading to the self-trapping of the synchrotron radiation, resulting 
in the growth of the radiation intensity, which manifests itself observationally as 
a flare. The condition for the self-trapping onset is derived, showing that if the 
phenomena takes place at lOOMeU, then it doesn’t happen at lower (or higher) 
energies. This specific kind of higher/lower energy cutoff could explain why when 
we observe the flare at 100 MeV no enhanced emission is observed at lower/higher 
energies! 

Subject headings: pulsars: general - pulsars: individual(PSR0531) 
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1. Introduction 

Crab nebula is a source of almost steady high-energy emission. Observations made 
from orbital probes (Fermi, SWIFT, RXTE) found an evidence of its variability in the X-ray 
range. Recently, Fermi and AGILE satellites detected brief and strong bursts of gamma 
radiation above 100 MeV, with its source located in the Crab nebula. One of such events 
was detected in September, 2010. Measurements made by AGILE on 12th and 16th of April, 
2011 have discovered new, unusually powerful bursts exceeding “September” ones more 


tha n five times. Finally, in 20 13, the most powerful bursts were observed flAliu et al. 2014J) . 


In (Buhler & Blandford 2014) an up-to-date review was given of the high-energy emission 
of Crab Nebula and the pulsar PSR 0532. 

The one and only source of radiation energy in the nebula is the rotational energy 
damping of the pulsar PSR 0532: 


A = m A = 5 ■ io 38 ^ 

at at sec 


( 1 ) 


where fl = 200sec~ 1 is the angular frequency of the pulsar and / = 10 45 g • cm 2 is its 
moment of inertia. This is more than enough to explain the value of the total luminosity 
of the nebula dW/dt = 5 • 10 38 erg/sec. Note that the luminosity of the pulsar radiation 
amounts only for about one percent of the total radiation intensity. In order to generate 


radiation with the energy up to 100 MeV, i 


is necessary to have particles wit h the energy 


of the order of PeV’s (10 15 eV). Tademaru ( Tademaru 19731) and W.J. Cocke f Cocke 1973h 


noticed that particles injected by the pulsar in the nebula, due to the radiation losses, 
can not have energies exceeding TeV range (10 12 ). The possible mechanism of the particle 
acceleration up to TeV ener gies and further production of 100 MeV photons was suggested 
in ([Clausen-Brown & Lyutikov 2012), where gamma-bursts were explained in terms of the 
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Doppler enhancement due to the variation of the averaged mini jet orientations in space. The 
latter, in its turn, happens due to th e magnetic field line reconnection in the region occupi ed 


by magnetically dominated p 


asrna (Blackman & Field 1994; 


Lyutikov & Uzdenskv 20031 


It was also shown (Kirk 2004) that in current layers, where E > B, while sufficiently fast 


particles move a 


ong the electric field lines, particles are e ntrapped by the layer and are 


accelerated. In (Cerutti et al. 2012a]; 


Cerutti et al. 2012bl ) a region was considered where 


the electric field exceeds the normal component of magnetic field E > B. It was shown 
that the mechanism of particle acceleration is similar to reconnection and for large-scale 
processes of the order of 10 16 cm. particles are accelerated up to PeV energies, which 
corresponds to Lorentz factors of the order of 3 x 10 9 . An obvious test of the proposed 
mechanism could be a considerable change of the structure of the magnetic field. So far 
such a noticeable change of this kind was not observed; in particular, Chandra failed to find 
such a variation in the X-ray range. 

Let us consider the possibility of the energy transfer to the nebula by the 
kinetic energy of particles. The mechanism of appearance of the particle acceleration 


( Klepikov 1954; 


Tademaru 1973 


Deutch 1955 

Erber 1966 

Ruderman 


dichel 1982)) 


Goldreich fe Julian 1969 


Sturrock 1971 


density of the electron flux is equal to mc 2 ^nGj, where Uqj - is the number density of 
particles ext racted from the pulsar surface by the electric hied, generated by the neutron 
star rotation ( Goldreich & Julian 19691 ). Near the light cylinder ugj = 10 6 cm~ 3 . While 


is the Lorentz factor of primordial particles. If we multiply this by the volume of the light 
cylinder and divide by the time interval, necessary for the transfer of the energy Rlc/ c i we 
find out that the power provided by the particle flux is equal to: 


rnc 3 nGj^bRic ~ 10 33 erg/sec. 


( 2 ) 
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Obviously, it is not sufficient. Therefore, it is necessary to have an additional 
mechanism of particle acceleration. This mechanism will be proposed in the next section 
of the paper. In the third section long wave-length electrostatic waves are suggested and 
considered as candidates for the energy transport, transferring energy from the pulsar to 
its nebula practically without any losses. In the fourth section we consider the process 


of part icl e acce 


eration in the Crab nebula related to the collapse of Langmuir waves 


(jZakharov 1972l) . The L angmuir wave collapse pheno menon was widely discussed in 


the plasma theory (e.g., flArcimowich fe Sagdeev 19791 )). but it is relatively unknown to 
astrophysical community. In order to explain observed lOOMeR gamma bursts in the fifth 
section we argue that their appearance is related to yet another plasma nonlinear process 
- self-trapping - which could explain unusual properties of these bursts. In the final, sixth, 
section, we summarize our results. 


2. Parametric generation of Lengmuir waves in the pulsar magnetosphere 


Recently it was shown flMachabeli et al. 20051: iMahaian et al. 201311 that relativistic 
centrifugal fo rce c an generate Langmuir waves. According to Machabcli & Rogava 


(iMachabel i & Rog ava 1994 ) when a particle moves along a straight, rotating field line the 


direction of centrifugal acceleration changes when particle’s initial velocity Vq < c/y/2 and 
the particle decelerates. It was shown that this regime of the motion leads to the parametric 
growth of Langmuir oscillations. 

Let us write down the equation of motion in an inertial frame of reference: 

+ (Vt ■ = -c 2 7 i(l - fl 2 r 2 ) 1/2 V(l - fl 2 r 2 ) 1/2 + q (e + ~^v x b\ (3) 


where p, E , B are momentum, electric held and magnetic held, respectively, is the 
angular velocity of rotation, r - distance from the center of the pulsar to the particle, 
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i = e~,e + - electron and/or positron. In the weak turbulence approximation, if we consider 
instantaneous values of p, E , B as sums of their regular and perturbational components - 
p = po + p, E = E 0 + E', B = B 0 - then we find that: 

dv 0 


= fl z r 1 


2 v o 

dt \ c 2 (l — fl 2 r 2 /c 2 ) 

which yields the following simple, periodic solution for relativistic velocities: 

Vq = c ■ cos (fit + 0) 

where 0 is an initial phase and c is the speed of light. 


( 4 ) 


( 5 ) 


We can write also the continuity equation and the Poisson equation, taking into 
account that n is the fluctuation of number density - deviation of this quantity from its 
equilibrium value no, v Q = pofm and v = p/m. This way we will obtain the closed set of 
equations for electrons and positrons. Further, let us consider a model of two streams with 
different Lorenz factors. One of them we denote with the index “p” and another with the 
index “b”. These two streams are connected with each other by means of the common 
electric field E '. Taking this fact into account and applying Fourie r tran sform we reduce 


the system of equations for both streams to the following equation (l Mah ajan et al. 2013): 

- 2 — r 


(<u 2 - (U 2 ) 


ryS 

Ip 


N p ( u) = — V' Ii(a)l s (a)exp[i(l(j) p -S(j) h )\ 
7 b ^ Ls 


(OJ 


N p (u-sfl+lfl ) (6) 


(c o — sh2) 2 

where N p (t) = n p exp[i ■ asin(flt + 0 P )], a = kc/fl, while I s (a) - are Bessel functions, 
s = 1, 2, 3,... and the sum is taken from — oo to +oo. Taking the average of (6) by phases we 
see that on the right hand side of the equation terms with 0 p = 0& give nonzero contribution 
when the harmonics coincide l = s. The equality of phases means that particles at the 
initial moment of time t = 0 are situated at the equal distance from the center of pulsar. 
The ultimate dispersion relation has the form: 


2 U P 

7 p 


cvz 


7 b 




or 


(c o — shl) s 


( 7 ) 
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Taking u = u 0 + A we find out that cOq = 2 cu 2 / 7 3 and ujq = sfl. The (7) relation contains 
two resonances which has to be satisfied simultaneously. For A we obtain cubic equation, 
with its solution containing an imaginary part: 

ImA = ±v / 3/2(2cn 0 a;^/f (a)/ 2) 1/3 ( 8 ) 


where ca 2 = 47 t e 2 Jib/m, the Bessel function index s = ca 0 /f2, ImA is the growth rate - 
the quantity which specifies the growth rate of electrostatic perturbation. This quantity 
depends on the quantitative value of the Bessel function. Its argument a = kc/fl » 1. For 
large values of the argument the Besse l function has m aximum value when the argument 
and index are equal s = a (Abramowitz & Stegun 1964). 


Therefore, when uq = kc , in the parametric interaction of particles with waves 
participate only very high-order harmonics with s » 1. It is important to know what 
fraction of the Langmuir wave energy is dissipated in the pulsar magnetospheric plasma. 
An effective energy dissipation mechanism is Landau damping, which takes place when 
v p h < v. For this purpose we must know at which wave numbers k c the dispersion curve 
crosses the c o = kc line. 


Electrostatic waves are described by the following equation: 

f dp 


i-E 


U P I 0 


(9) 


(cu — kv) 

Here the summation is made by the particle species - plasma electrons and positrons, 
ujp = 87 Te 2 n p /m, while / is the plasma distribution function normalized on unity. The (9) 
equation can be easily solved in two extreme cases: c o » kv and ca ~ kc. The first case 
describes Langmuir waves, with phase velocities is exceeding the speed of light and 

cu 2 = w 2 7 -3 + 3 k 2 c 2 (10) 

For the second case in (ca — kv) 2 we substitute co ~ kc and v = c(l — 7 -2 ) -1 / 2 . If we expand 
the quadratic root by the small parameter I/ 7 2 from ( 10 ) we derive the expression for the 




value of wave number vector k c , at which the curve crosses the u = kc line. 


co 2 /c 2 ee kl = 2'yuj 2 /c 2 


and the spectrum has the followng form flLominadze et ah 19791 ) 


co = kc[ 1 — a{k — k c )/k\ 


(ii) 


( 12 ) 


From the ( 8 ) condition it is evident that co 0 < kc, the co 0 = kc condition is fulfilled only 
when k < k c . Therefore, the growth rate is maximum in the region where the phase velocity 
ooo/k > co 0 /k c ~ c. But there are no particles with such velocities, hence Langmuir waves 
can interact with the particles of electron-positron plasma only via nonlinear effects. 


We have to note that while the parametric instability is developed the energy of the 
pulsar rotation is directly pumped into the energy of Langmuir waves. From (4) it is easy 
to estimate potential power of the parametric instability: 


dW 

dt 


co 2 r 3 c 3 (m e n LC ) 
(c 2 — fl 2 r 2 ) 


erg/sec 


(13) 


While the longitudinal velocity v\\ is decreased the linear velocity of rotation fir is increased 
and consequently W t may reach the value 5 • 10 38 at the light cylinder. Hence, we see 
that due to the parametric interaction, strongly relativistic particles manage to pump the 
rotational energy of the pulsar directly into the Langmuir waves. 


Naturally, the question arises as to what causes the upper limit to the growth of 


the luminosity with the value 5 • 10 3S erg 

reconstruction of the magnetic field from c 

/sec. In our opinion the limit is relat 

ipole to the monopole configuration ( 

ed to the 

Michel 1982; 

Rogava et al. 2003: 

Osmanov et ah 2008; 

Osmanov et ah 2009: 

Bogovalov 200; 

D- 


We found out that the parametric generation of Langmuir waves for c o/k > c happens 
with higher growth rate than for a i/k < c. Therefore, Langmuir waves do not damp. On 
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the contrary, the Langmuir oscillations become accumulated and they transfer energy along 
the magnetic held lines up to the Crab nebula. This process is considered in detail in the 
next chapter. 


3. Transfer of Langmuir oscillations to the nebula 


Let us make sure that the energy losses are negligible when the waves are propagating 
towards the Nebula. It was assumed that the most probable mechanism of longitudinal 
wave-particle interaction, in the weak turbulence approximation, is a wave scattering by 
plasma particles. This statement is true if the scattering process is due to the resonance of 
the type of 

{u k ~uj kl )/(k - h) = v (14) 

ffere, if we take values of the frequencies from (10) and (12), then it appears that the waves 
with the spectrum cu ~ fee are scattered onto the thermal particles of relativistic plasma, 
causing the energy pumping into the long-wave band. 

The plasma oscillations can be considered as a gas of quantum quasi-particles - 
plasmons - (analog of photons) and the concept of the number of plasmons is introduced. 
Scattering does not change the number of plasmons, which is a conserved quantity: 

N k = \E k \ 2 /uJk — const (15) 


The energy is pumped to the long-wavelength region of k —> 0. From the conservation of 
\E k \ 2 /uJk it follows that \Ek=o\ 2 /\E k \ 2 ~ u>k=o/^>k ~ 7~ 2 - Thus, th e wav e energy losses are 
quite significant. This result has been obtained by (iLom i nadze et ah 1 9791 ). based on the 
assumption that the waves are generated in the subluminal region, where v p h < c. 


ffowever, in our case, the main bulk of energy is generated in the region where 
Therefore, the nonlinear scattering of Langmuir waves onto the plasma particles should 



10 


be considered in the region where v p h S> c. Substituting the value of electrostatic wave 
frequency from equation (10) in the resonant condition of nonlinear scattering, we obtain 
the following expression: 


_ uj p /l + 3fc 2 c 2 7 3 1 + 3/c' 2 c 2 7 3 \ 1 

7 p 3 / 2 \ 2a j p 2 2a j p 2 ) k — k' 

Taking into account that k 2 — k' 2 — (k — k'){k + k') and assuming k ~ k ', we 
derive the following estimation: v/c ~ 3c/v p h 1. The thermal particles of this sort in 
ultra-relativistic plasma practically do not exist and therefore the scattering is not effective. 

Since the density of the plasma particles is decreasing as the wave propagates towards 
the Nebula, co p 2 decreases as well. However, co 2 should be conserved. In order for this 
to happen, in the non-relativistic case, it is necessary the raise the value of k 2 . But in 
the relativistic case, the frequency approximately equals cu p ~ oJlc/^lc 3 ^ 2 when k —> 0. 
Therefore, it is sufficient to have: 



1/2 



1 


(17) 


Here, ulc ~ 10 13 cm” 3 is the bulk plasma density at the pulsar light cylin- 
der; n-w ~ 10 3 cm” 3 deno tes the plasma density within the fibers of Crab Nebula 
({M anchester & Taylor 19771) . The frequency does not change while Langmuir waves are 
propagating towards the nebula. Therefore, as it follows from the condition (17), the 
thermal particles of plasma have to be cooled down to Lorentz factors of order ~ 1 — 10. 

It means that, due to the ineffective wave-particle interaction, the energy of the long waves 
(k —> 0) is transferred to the distance, Rn ~ 10 17 cm, almost without any considerable 
losses. 


In the Nebula, the spectrum of the Langmuir oscillations is also situated beyond the 
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light cylinder k —> 0. The energy of electrostatic oscillations has to be converted into the 
kinetic energy of the particles. This issue is considered in the following section. 


4. Collapse of Langmuir oscillations 


“Condensate” -is an intense gas of long-wave plasmons which is unstable when the 
background is modulated by low-frequency waves. But there are no such waves in the 
pulsar magnetosphere, since the magnetosphere consists exclusively of electrons and 
positrons. In the Crab nebula, unlike the magnetosphere, there are ions and they can 


ensure the low-frequency modulation of the p lasmon 


background, which, in its turn, leads 


to modulational instability ( Vedenov & R.udakov 1962 ). The development of this instability 


manifests itself by the formation of caverns - regions of localization of excessive wave energy. 
Due to this instability the energy of Langmuir waves is localized in the cavern. Plasmons 
appear to be locked in the cavern. High-frequency pressure pushes particles out of the 
cavern. Consequently the cavern deepens and it absorbs in itself more and more plasmons. 
The process of spatial attraction of plasmons is accelerated. It is accompanied by the 
collapse of the cavern. 


This phenomenon can not be considered in the framework o 


the weak turbulence 


approximation. The dynamics of cavern collapse was studied in (iGaleev et al. 19771) . A 
numerical ex pe riment, illustrating the dynamics of the cavern collapse was presented in 


flDegtiarev et al. 19761 ). It turned out that in the process of cavern collapse the wave vector 


k grows until the collapse phase velocity equals the speed of particles co/k = v. Particles 
which fall in the resonance with the wave extract the wave energy. Due to Landau damping 
the cavern colla pse stop s to develop and the distribution function of the particles acquire a 
long tail ( Pelletier 1982 ). 
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The above-described collapse scenario takes place in the plasma without external 
magnetic field. In Crab Nebula plasma number density is no = 10 3 and magnetic field 
B 0 = 1CT 3 G, implying that co p » u>b = eB 0 /mc and the influence of the magnetic field 
on the collapse of Langmuir waves can be neglected. In this case the equation for the 
amplitude of Langmuir oscillations can be written in the following way: 


d 2 E 

dt 2 


“'p_ 

3/2 

7 p 


E — 


3 rf,u‘ d 2 E uj- 


r jp/ 2 dx 2 


3/2 

7 p 


n , 

— 1 E 

n 0 


(18) 


At the other hand the number density perturbation, resulting from the plasmon pressure, 
obeys the equation: 


where M is the ion mass. Auto-model solution of the system have the “burst-like’' 
appearance: 


(19) 


\E\ 2 = 


El 


(To - t ) 2 


( 20 ) 


where r = n 0 T e (uj p /7o / 2 )tEl/Sit. 

Here | _£7 0 1 2 is the held energy density at the initial 
and r 0 - is the time of the active phase of the collapse flZakharov 1972( 1. In the course of the 


stage o f the modulational instability 


collapse, when waves are being damped, electric field is decreasing with the increase of the 
particle kinetic energy. 


5. Synchrotron radiation of the nebula and 100 MeV bursts 


As early as in 1953 I.S. Shklovsky in order to explain the continuous radio emission of 


the Crab Nebula suggested the synchrotron mechanism ( Shklovslv 19531 ). Soon afterwards 


observations in the visible spectrum range revealed a substant 

(Vashakidze 1954; 

Dombrovskv 1954: 

Oort & Warvalen 1956) 


ial level of linear polarization 
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polarization degree in various regions evidently indicates that synchrotron mechanism has 
indeed the leading role in the generation of the nebula radiation. 


In ( Gould 1965 


Weekes et ah 1989 


de Jader et ah 1 9961 ) the mechanism of Compton- 


synchrotron formation of spectra was proposed. Above a few hundreds of MeV the intensity 
of synchrotron emission is strongly diminished and starting from frequencies of the order of 
1 GeV the radiation is determined by the inverse Compton scattering. In this paper we 
will restrict ourselves by considering frequencies of the order of lOOMev, which means we 
consider only the synchrotron mechanism of the radiation. However, we have to note that 
the problem of the pulsar rotational energy transfer to the Nebula is open and actual also 
for the inverse Compton radiation. Not rejecting the idea of explanation of gamma bursts in 
terms of the magnetic field line reconnection, which was mentioned in the Introduction, we 


Chiso et ah 1961: 


on the medium (i 

\skarvan 1962: 

'l -- | 

Whitham 1974; 

Akhmanov et al. 1967: 

Klimontovich 196 

3). As it turns out, bursts with lOOMeV energies 


are not accompanied by any tangible in crease of radiation intensity in any other ranges 


of the Crab Nebula radiation (lAliu et al. 2014). Such behavior of the bursts makes us 


surmise that this phenomenon can be explained by the self-trapping effect. In this case the 
mysterious behavior of the bursts could be explained in a quite natural way. 


Self-trapping is a well studied phenomenon in nonlinear optics, its theory was 
developed in 1960s, accompanying the appearance of powerful light sources - lasers. As it 
was mentioned above in Crab Nebula a modulational instability is developed which leads 
to the appearance of caverns, where the energy of long wavelength electrostatic oscillations 
is stored. Caverns are eventually collapsed. At the beginning of the process the collapse 
time scale is determined by the growth rate of the modulational instability. Further, while 
r —» tq, according to (20), the collapse rate increases in a burst-like way; the active phase of 
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the collapse starts and Landau damping makes substantial the absorbtion of the energy by 
cavern particles from Langmuir waves. Particles get accelerated up to required limits and 
the cavern collapse stops. The size_of the cave rn throughout its active phase is estimated as 


a few Debye radii wide (lArcimowich fe Sagdeev 19791 1. As it is well-known, Debye radius is 


the scale of species separation in plasma and is determined by the length scale, over which 
electron density perturbation may be shifted due to its thermal motion over the period of 
plasma oscillations. In our, relativistic, case: 


r D = 


( 21 ) 


For the bulk of the plasm a par ticles in th e Crab Nebula the particle number density is of 
the order of n 0 10 3 ([Manchester fe Tavlor 19771 ). Prior to acceleration Lor entz factors 
of plasma particles are of the order of y p « 1 — 10 ( Machabeli fc Usov 19891 ). After the 


acceleration happens, the energy maximum, as it follows from observations, falls in the 
soft X-ray range. That is why one can assume that the majority of plasma particles are 
accelerated up to Lorentz factors 7 ^ 5 - 10 5 . Then the frequency of synchrotron radiation 
is situated in the right range: 


cu s = (eB/mc) y 2 = 1.8KeV 


( 22 ) 


Putting the values of the number density and Lorentz factor in (21) we obtain that 
ro — 5 x 10 14 cm. An average distance between plasma particles is r cp n -1 / 3 = 10 _1 , hence 
we can estimate maximum quantity of particles in Debye volume as Np ^ {rr>/r cp ) z 10 46 . 

In the process of the cavern collapse participate also particles from the distribution 
function tail, which are being accelerated up to Lorentz factors 7 « 3 x 10 9 . From (22) we 
can estimate that particles with such Lorentz factors are emitting photons with energies 
about lOOMeD due to the synchrotron radiation mechanism. Assuming that n p y p = n s ^ sy 
we can find the number of particles N s ~ 10 38 . It is close to the value found in observations 
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flAliu et ah 2014). Multiplying the number of particles on the radiation intensity for one 
particle: 


de 

dt 


2e 2 u> B 

3c 


7 


(23) 


we obtain the intensity of the synchrotron radiation of the order of lOOMeV. 

g|) « 10 x erg/a 


(24) 


The electromagnetic wavelength can not be less than r B since in the Debye volume 
electromagnetic field is screened out due to the grouping of charged particles, that is, the 
polarization of the medium takes place. For convenience, let us imagine the following model: 
the medium within the Debye radius consists of dipoles with the same spatial orientation. 
In this case the polarization vector: 

P = N(ef) (25) 

In the course of the collapse transverse size of the cavern (in the X and Y directions) 
is decreasing. The acceleration of electrons happens in the same directions, which, in turn, 
generates synchrotron photons. Let us assume that local magnetic field is directed along 
the Z axis and the observer is located in the X direction. Synchrotron radiation is directed 
in the X and Y direction. The electric field of the radiation E(t) shifts charged particles. 
The shifting causes the appearance of the elastic counterforce /(f) = —r)f(t), where ij is the 
elasticity coefficient. For not too small values of E(t ) the elasticity force has nonlinear form: 


/(f) = -Jff(i) - qf(ty 


(26) 


The value of the electron shifting is determined from the equation of motion. In the 
relativistic case it has the following form: 

nd 2 r dr _ „ 

m'y —- — mV ——b rjr + qr = eE 
dt 2 dt 


( 27 ) 
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where the second term on the left hand side defines the dissipation - damping rate. Taking 
into account the definition of the polarization vector (Eq. (25)) and r) = m.'yujQ we obtain: 

+ T^y- + (cUo/ 7 2 )-P+ ( q/me 2 N 2 '-y 3 )P 3 = e 2 NE/m q 3 (28) 

In our case oj q is the frequency of Langmuir oscillations u a = ca^/q 3 / 2 . The E held is 
large but still it is much less than the internal held of the cavern. In particular, the power 
of the wave energy of the cavern stored in Langmuir waves is of the order of 10 38 erg/s, 
whereas the synchrotron radiation intensity is of the order of 10 36 erg/s. Consequently, the 
electric held of the wave is small compared to the internal held. In such circumstances the 
nonlinear term can be considered small and the equation can be solved by the method of 
consecutive approximations. Splitting P = Pl + Pnl, with Pl » Pnl and neglecting the 
nonlinear term we get: 

+ + ( w o/7 2 )^ = e 2 NE/m q 3 (29) 

Taking E = Acos(ut) we can find the solution: 

P{t) = e 2 NAcos(cut + <f>)/mq 3 ^/(w 2 - u$) 2 + 4T 2 o ; 2 (30) 

where tg$ = Tu;/( u 2 — ujq). 

Note that the Langmuir frequency is much smaller than the synchrotron radiation 
frequency ojq « co. Therefore, the range of frequencies which we are considering is far 
from the resonance: \uJq — u 2 \ » 4T 2 . It means that the contribution related to the 
dissipation can be neglected. The polarization vector P is related to the electric held via 
the polarizability of the medium /i, in the following way: P = / iE . This expression can be 
written also as: 


P(t) = n(uj)E{t) 


( 31 ) 



-17- 


The nonlinear approximation equation has the following form: 

' j t2 + UqPnl = qn 3 (u)E(t) /m^e 2 N 2 

(32) 

Let us rewrite E 3 (t) using the trigonometric identity: cos 3 (cut) = (l/A)(cosujt + cos3ut). 

This way on the right hand side of (32) we have two terms describing input of the first and 

the third harmonics. Consequently we can write: 

P{t) = fx(u,A)E(t) 

(33) 

n(cu,A) = /x(u;)[l + 3q/d 2 (co)A 2 /4mn 2 e 2 uj 2 ] 

(34) 


and h(uj) is determined from the solution of (29): 

n(u) = e 2 N/m r y 3 io 2 (35) 

Note that in the series expansion of A) we retain first non-vanishing terms. 

The dielectric permittivity of the medium is described by the tensor £ij(uj,E). The 
connection between this tensor and /j,ij(uj,A) is determined by the expression: 

£ij(uj, E) = Sij + 47 A) (36) 

The induction vector D — E + 47rP, P,; = £ij(u, E)Ej. After this is written down, 
taking into account (36) and (34), we write down Maxwell equation: 

_ f)E 

V x B = (1 /c)[e{(jS) + 37Tg/i 3 (w)4 2 /mAf 2 e 2 w 2 ]— (37) 

From this analysis it is clear that the influence of the nonlinear term is equivalent of 
the change of the dielectric permittivity or the refraction index of the medium. When an 
electromagnetic wave is propagating in the medium the refraction index H = c/v p h is a 
function of frequency. Therefore, the dispersion of light depends on the refraction index. 
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From the expression H 2 = e we find that in the cavern, which undergoes the active phase, 
the refraction index becomes equal to H = H L + H NL , where 

H 2 L = e(u), H nl = H 2 A 2 , H -2 = 6nqii 3 (u)/mN 2 e 2 uj 2 (38) 

Therefore, if H 2 > 0, the refraction index in the cavern H = Hl + Hnl turns out to be 
larger than the refraction index beyond the cavern, which is equal to H = Hl- In the whole 
volume of the cavern let us separate rays which are directed towards the observer. These 
rays are generated by particles with Lorentz factors 7 = 3 x 10 9 and directed towards the 
observer. That’s why they are emitted within the narrow cone with the opening angle I/ 7 . 
Due to the linear diffraction they have to diverge - diffuse in the direction normal to the 
propagation direction and before leaving the cavern they will be located within the cone 
with opening angle 2 9d, where 9d ~ X/tdHl- However, when these rays leave the nonlinear 
medium of the cavern and enter the space with refraction index Hl they undergo nonlinear 
refraction. If a ray falls on the boundary between nonlinear, optically thicker medium and 
linear, optically thinner medium and if the incident angle satisfies the condition 9q > 9d 
then all diffracted rays will undergo complete internal reflection. We are interested in the 
condition 9 0 & 0 D which leads to the situation when the rays are combined to form a 
parallel beam and the observer witnesses an increase of the radiation intensity. 

An asymptotic sliding angle for the complete reflection from the cavern boundary is 
determined in the following way: cos9$ = Hl/(Hl + H 2 A 2 ). For the small 9 angle we find: 

9 0 « 2 (H 2 /H l )A 2 (39) 

From the condition 9q ~ 9d we determine the length of those waves which are focused in a 
parallel beam: 

A « 2A 2 H 2 r D (40) 

Substituting (35) and (38) we find out that H 2 ~ 1/ca 8 and if the self-trapping 
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condition is satisfied for frequencies close to lOOMeV, from the o bserved c 


radiation intensity on freq uency (see, e.g., Fig. 1 in fide Jader et al. 19961) ) or Fig. 2 in 


ependence of the 


(Buhler & Blandford 2014)) we conclude that for other frequencies the (37) condition is not 
met. This circumstance could explain strange behavior of bursts in the 100 MeV region of 
the spectrum, while at the same time in other frequency ranges the increase of the intensity 
is not observed. 


6. Conclusion 

In the conclusion we just briefly summarize the contents and main results of this paper: 

• The possibility of the direct pumping of energy from the neutron star’s vast rotational 
kinetic energy storage 5 • 10 38 erg/sec to proper electrostatic plasma (Langmuir) 
oscillations is demonstrated. 

• It is shown that the growth rate of the perturbations is maximum in the “superluminal” 
area, where phase velocity of perturbations is exceeding the speed of light. That is 
why in this region the condensate of plasmons is formed, which is transferred to the 
Crab Nebula. 

• It is demonstrated that the transfer of the energy of the plasmon condensate from 
the pulsar magnetosphere to the nebula over the huge distance 3 • 10 17 cm takes place 
practically without any tangibles losses. 

• Unlike the pulsar magnetosphere in the nebula apart from electrons and positrons 
there are also protons. That is why a modulational instability is developed, which 
leads to the collapse of Langmuir waves. A cavern is formed which collapses and on 
the final stage of the collapse particles attain very high Lorentz factors, resulting in 
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the powerful emission of the nebula and the collapse stops at the scale of a few Debye 

radii. 

• Theoretical estimation of the cavern (radiating region) and the number of emitting 
particl es co incides with evaluations made on the basis of observational results in 
(Ahu ct ah 2014 b 

• It is shown that in the course of the active phase of the collapse in the cavern, due 
to the influence of nonlinear processes on the polarization properties of the medium, 
self-trapping of the synchrotron radiation generated within the cavern, takes place. 

• It is shown that if the conditions for the appearance of self-trapping phenomenon are 
fulfilled for certain values of emitted wave frequencies, for other, both higher and 
lower values of the frequency, they are not satisfied. 
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